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ABSTRACT 

Probably No! As an example, using soft EOSs consistent with existing terrestrial nuclear laboratory 
experiments for hybrid neutron stars containing a quark core described with MIT bag model using 
reasonable parameters, we show that the recently discovered new holder of neutron star maximum 
mass PSR J1614-2230 of 1.97 ± O.O4M0 can be well described by incorporating a Yukawa gravita- 
tional correction that is consistent with existing constraints from neutron-proton and neutron-lead 
scatterings as well as the spectroscopy of antiproton atoms. 
Subject headings: hybrid star, hyperon, quark, gravity 



1. INTRODUCTION 

What is gravity? Are there additional spacetime di- 
mensions? These are among the Eleven Science Ques- 
tions for the New Century identified by the Commit- 
tee on the Phy sics of the Univ erse, US National Re- 
search Council (|CPUNRar2003| ). Interestingly, despite 
the fact that gravity is the first force discovered in 
nature, the quest to unify it with other fundamental 
forces remains elusive because of its apparent weakness 

(lArkani-Hame d et al 



at short-distaiice, see, e.g. , refs. ^^ ..^.j.... 
199 1 |Pea|[200ltlHov leli2003l: 'Lo ng et al..,200a:Jean et al 
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veloping grand unification theories, the conventional 
inverse-square-law (ISL) of Newtonian gravitational force 
has to be modified due to either the geometrical ef- 
fect of the extra spacetime dimensions predicted by 
string theories and/or the exchange of weakly interact- 
ing bosons, such as the neutral spin-1 vector /7-boson 
(|Fayct 2009,) , proposed in the super-symmetric extension 
of th e Standard Model, see, e.g., refs . (Adclbcrgc r et al.l 
20031 [20091 iFischbachI 11991 INewman 2009: Uzani l2003r 
RevnaudI 12001 " for recent reviews. The modified 
gravity has also been proposed as an explanation for 
the presen t period of cosmq logical acceleration, see, 
e.g., ref. pePeo et al.l 120081 ). The search for evi- 
dence of modified gravity is at the forefront of re- 
search in several sub-fields of natural sciences includ- 
ing geophysics, nuclear and particle physics, as well 
as astrophysics and cosmo l ogy, see, e.g., refs. 
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Azam et al.l 120081: iGeraci et al.l ~ l2010t iLucchesi et al l 
20101) . Various upper limits on the deviation from the 



ISL has been put forward down to femtometer range. 
Since the composition of neutron stars are determined 
mainly by the weak and electromagnetic forces through 
the P equilibrium and charge neutrality conditions while 
their stability is maintained by the balance of strong and 
gravitational forces, neutron stars are thus a natural test- 
ing ground of grand unification theories of fundamental 



forces. Moreover, neutron stars are among the dens- 
est objects with the strongest gravity in the Universe, 
making them ideal places to test strong-field predictions 
of General Relativity (GR) (Psahis 2008). The masses 
and radii of neutron stars are solely determined by both 
the strong-field behavior of gravity and the Equation of 
State (EOS) of dense stellar matter. However, there is no 
fundamental reason to choose Einstein's equations over 
other alternatives and it is known that the GR theory 
itself may break down at the limit of very strong gravita- 
tional fields, see, e.g., ref. (Psaltis 2008) for a comprehen- 
sive review. In fact, effects of modified gravity on prop- 
erties of neutron stars have been under intense investi- 
gation. As expected, results of thes e studies are strongl 
model dependent, see, e.g., refs. (Germani et al.' '2001J 
Wiseman 2002; Azani et al.l [2008; Krivoruchcnko et al 
20091: IWen et al.li2009D . Nevertheless, it is very interest 



ing to note that alternative gravity theories that have 
all passed low-field tests but diverge from GR in the 
strong-field regime predict neutron stars with signifi- 
cantly different pro perties than their GR counterparts 
(jPePeo et al.|[2003t ). Moreover, the deviations for neu- 
tron star properties from the GR predictions for these 
theories are larger than the uncertainty due to the poorly 
known EOS of dense matter in neutron stars. It was also 
clearly shown that the neutron star maximum mass alon e 
can not distinguish gravity theories (jPePeo et al.ll2003l) . 
Furthermore, in the endeavor of testing GR theory of 
gravity using properties of neutron stars, it is known 
that there is a degeneracy between the matter content 
and gravity. This degeneracy is tied to the fundamental 
Strong Equivalence Principle and can only be bro ken by 
using at least two independent observables ( Yune s~et al.l 
12010) . 

Recently, using the general relativistic Shapiro delay 
the mass of PSR J 1614-2230 was p recisely measured to 
be 1.97 ± 0.04M,7. (|Pemorest et al.fl2010[) . making it the 
new holder of the maximum mass of neutron stars. Com- 
paring with mass-radius relations predicted from solving 
the TOV equation using various EOSs within GR theory 
of gravity, it was shown that the mass of PSR J1614- 
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2230 can rule out almost all soft EOSs especially those 
associated with hyperon or boson condensation. While 
conventional quark stars with soft EOSs are also ruled 
out by this observation, neutron s tars with strongly in- 
teracting quark cores are allowed (jPemorest et al.ll201Cll : 
lOzel et alJ 120101: iLai et al.|[20Tol ). It was further shown 
that a transition to quark matter in neutron star cores 
can occur at densities comparable to the nuclear satu- 
ration density po only if the quarks are strongly inter- 
acting and are color superconducting (|Ozel et al.ll2010D . 
The mass of PSR J1614-2230 was then used to constrain 
the interacting parameters of quarks. It was also shown 
that neutron stars with interacting quark clusters in their 
cores or solid quark stars can be very massive. Using the 
Lennard- Jones potential for interactions between quark 
clusters, the mass of the PSR J1614-2230 was used to 
constrai n the nurnber o f quarks inside individual quark 
clusters (i Lai et al]|2010f) . In this work, using soft nuclear 
EOSs for hybrid stars containing a quark core described 
by the MIT bag model with reasonable parameters, we 
show that the mass of PSR J1614-2230 is readily ob- 
tained by incorporating a Yukawa gravitational correc- 
tion that is consistent with existing constraints from ter- 
restrial nuclear laboratory experiments. 

2. NON-NEWTONIAN GRAVITY AND MODEL EOS OF 
HYBRID STARS 



Fujii (jFuiiil |1971|) first proposed that the non- 
Newtonian gravity can be described by adding a Yukawa 
term to the conventional gravitational potential between 
two objects of mass mi and m2, i.e.. 



Gmirn2 
V{r) = (1 
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where a is a dimensionless strength parameter, A is the 
length scale and G is the gravitational constant. In 
the boson exchange picture, a = ±g^/(47rG'm^) where 
± stands for scalar/vector bosons and A = l//i (in 
natural units). The and /i are the boson-baryon 
coupling constant and the boson mass, respectively. 
The light and weakly interacting J7-boson is a favorite 
candidate mediating the extra interac tion (jFaveti 120091 : 
iKrivoruchenko et al. J [20091: [Zh{j|2007[ ). Similar to the 
degeneracy between matter content and gravity, there 
appears to be a duality of incorporating effects of the 
Yukawa term in either the TOV equati on or the input 
EOS. Nevertheless, according to Fujii (Fujii' '1988), the 
Yukawa term is simply part of the matter system in gen- 
eral relativity. Therefore, only the EOS is modified and 
the TOV equation remains the same. Within the mean- 
field approximati on, the extra energy density due to the 
Yukawa term is (jLong et al.l 120031 : IKrivoruchenko et al..l 
[200i) 
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where V is the normalization volume, p is the baryon 
number density and r = \xi — X2\- Assuming a con- 
stant boson mass independent of the density, one obtains 

the corresponding addition to the EOS Pub = ^^P^j 
which is just equal to the additional energy density. As 
it was emphasized by Fujii (Fuiii .1988,) . since the new 
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Fig. 1. — (Color online) Model EOSs for hybrid stars with 
and without the Yukawa contribution using MIT bag constant 
BV4=170 MeV and B^l^=\'m MeV, respectively. 



vector boson contributes to the EOS only through the 
combination g^/^"^, while both the coupling constant 
g and the mass /i of the light and weakly interacting 
bosons are small, the value of j p^ can be large. On 
the other hand, by comparing with the g^ j [i? value of 
the ordinary vector boson w, Krivoruchenko et al. have 
pointed out that as long as the value of the U 

boson is less than about 200 GeV"^ the internal struc- 
tures of both finite nuclei and neu tron stars will not 
change ([Krivoruchenko et al..l [2009| ). However, global 
properties of neutron stars can be significantly modified 
(Krivoruchenko ct al., 2009; Wen ct al. 2009) One of the 
key characteristics of the Yukawa correction is its compo- 
sition dependence, unlike Einstein's gravity. Thus, ide- 
ally one needs to use different coupling constants for var- 
ious baryons existing in neutron stars. Moreover, to our 
best knowledge, it is unknown if there is any and what 
might be the form and strength of the Yukawa term in 
the hadron-quark mixed phase and the pure quark phase. 
Nevertheless, instead of introducing more parameters, for 
the purpose of this exploratory study, we assume that the 
Pub term is an effective Yukawa correction existing in 
all phases with the g considered as an averaged coupling 
constant. 

Including the Yukawa term the EOS becomes P = 
^0 + Pub where Pq is the conventional pressure inside 
neutron stars. For the latter, we use typical model EOSs 
for hybrid stars containing a quark core covered by hy- 
perons and leptons. The quark matter is described by the 
MIT bag model with reasonable parameters widely use d 
in the literature ([Chodos et al.1 1 1974) : [Heinz et al.lll986( ). 
The hyperonic EOS is modelled by using an extended 
isospin- and momentum -dependent effective interaction 
(MDI ([Das et al.ll2003[ )) for the baryon octet with pa- 
rameters constrained by empirical properties of symmet- 
ric nuclear rnatter, hyper-atoms, and heavy-ion reactions 
([J. Xu et al.ll2010( ). In particular, the underlying EOS 
of symmetric nuclear matter is constrained by compar- 
ing transport model predictions with data on collective 
flow and kaon production in rclativistic heavy-ion colli- 
sions (Daniclewicz et al... 2000; .Li et al...2008X Moreover, 
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Fig. 2. — (Color online) The mass-radius relation of static neu- 
tron stars with i?^/'*=170 MeV and various values of in units 
of GeV~^ indicted using numbers above the lines. 

the nuclear symmetry energy Esym{p) with this inter- 
action is chosen to increase approximately linearly with 
density (i.e., the MDI interaction with a symmetry en- 
ergy parameter x=0 ()J. Xu et al.l 120100 ) in agreement 
with availabl e constraints aroun d and below the satura- 
tion density ()C. Xu et al.ll2010l ). It is worth noting that 
this Esym {p) is very simil ar to the well-known APR pre- 
dictio n up to about 5po ()Akmal et al.l 119981 : IXiao et ahl 
120091) . The Gibbs construction was adopted to describe 
the hadron-quark phase transition (Glcndcnning 200l]). 
Similar to the previous work in the literature, the hy- 
brid star is divided into the liquid core, inner crust and 
outer crust from the center to surface. For the inner 
crust, a par ameteriz ed EOS of Pq = a + he"^/^ is adopted 
as in refs. ( J. Xu et al. 2009). For the outer crust, the 
BPS EOS jSavm et al.lll971l) is adopted. As an exam- 
ple, shown in Fig. |l] are the EOSs for hybrid stars with 
MIT bag constant B'^/'^ = 170 MeV and B^'^ = 180 
MeV, respectively, with and without the Yukawa contri- 
bution. The corresponding hadron-quark mixed phase 
above/below the pure hadron/quark phase covers the 
density range of p/po = 1-31 to 6.56 and p/po = 2.19 to 
8.63, respectively. Including the Yukawa term the EOS 
stiffens as the value of / p^ increases. It is seen that the 
two sets of parameters, i?^/* = 170 MeV and g^/p^ = 
50 GeV-2 or B^'^ = 180 MeV and gV^^ = 40 GeV"^ 
lead to approximately the same total pressure. As a ref- 
erence, the MDI EOS for the npep matter is also shown. 

3. MAXIMUM MASS OF HYBRID STARS WITH 
NON-NEWTONIAN GRAVITY 

As an example, shown in Fig. [5] is the mass-radius 
relation of hybrid stars with the bag constant i3^/'*=170 
MeV and varying values of 5^/p^. First of all, without 
the Yukawa contribution (black solid line) the maximum 
stellar mass supported is only about 1.46Mq. Including 
the Yukawa term, as the EOSs are increasingly stiffened 
with larger values of 5^/p^, the maximum stellar mass 
increases. With = 50 GeV~^ the maximum mass 

of 1.97Mq is just in the middle of the measured mass 
band of PSR J1614-2230. The corresponding radius is 
about 12.4 km. To see more clearly relative effects of the 
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Fig. 3. — (Color online) The maximum mass of neutron stars 
as a function of g^/^a ^j^jj qI/a^ijq y^^y ^nd B^/^=im MeV, 
respectively. 
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Fig. 4. — (Color online) Constraints on the strength and range 
of the Yukawa term from terrestrial nuclear experiments in com- 
parison with Ri 40 — 50GeV~^. 

bag constant B and the Yukawa term, shown in Fig. [3] are 
the maximum stellar masses as a function of (7^/p^ with 
Bi/4 ^ MeV and B^/^ = 180 MeV, respectively. 
As expected, with B^^^ = 180 MeV a smaller value of 
g^/ p^ = 40 GeV~^ is needed to obtain a maximum mass 
consistent with the observed mass of PSR J1614-2230. 

4. COMPARISON WITH TERRESTRIAL CONSTRAINTS ON 
NON-NEWTONIAN GRAVITY AT SHORT DISTANCE 

As mentioned earlier, significant efforts have been de- 
voted to constrain the possible non-Newtonian grav- 
ity using terrestrial experiments. These experiments 
have established a clear trend of increased strength a 
at shorter length A. In the short range down to A « 
10~^^ — 10~^ m, neutron-proton and neutron- lead scat- 
tering data as well as the spectroscopy of antiproton 
atoms have been used to set upper limits on the value of 
g^/p^ (or equivalently the \a\ vs A). It is thus interesting 
to compare the values of 5^/p^ necessary to support the 
PSR J1614-2230 within the model presented above with 
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the constraints extracted from terrestrial experiments. 
While the range parameter A is expected to be much 
larger (smaller) than the radii of finite nuclei (neutron 
stars) , the maximum mass of neutron stars alone is not 
sufficient to set separate constraints on the values of a 
and A. Shown in Fig. |4]is a comparison with the terres- 
trial constraint s (Kamvshk ov et al.l 120081: [Barbieri et al.l 
[1971 lPokotiio^skit.200& : ,Nesvizhevskv et al.ll2008[ )" in the 
\a\ vs A plane. The straight line is for g'^/fi^ « 40 — 
50 GeV~^. ft is seen that the values of .g^//i^ neces- 
sary to describe the maximum mass of PSR J1614-2230 
are consistent with the upper limits from the terrestrial 
experiments. 



5. CONCLUSIONS 

Among all fundamental forces, gravity remains the 
most uncertain one despite being the first discovered in 
nature. Neutron stars are natural testing grounds of 
grand unification theories of fundamental forces. In par- 
ticular, they are ideal places to test GR predictions at the 
strong-field limit. Interpretations of observed properties 
of neutron stars require a comprehensive understanding 
of both gravity and the EOS of dense stellar matter. Be- 
fore strong-field gravity is well understood, it is unlikely 
that the maximum mass of neutron stars alone can rule 



out any EOS. As an example, using soft nuclear EOSs 
consistent with existing terrestrial experiments for hy- 
brid stars containing a quark core described by the MIT 
bag model with reasonable parameters, the maximum 
mass of PSR J1614-2230 is readily obtained by incorpo- 
rating the Yukawa gravitational correction that is con- 
sistent with existing constraints from terrestrial nuclear 
laboratory experiments. 
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